Edge States in Doped Antiferromagnetic Nano-Structures 
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We study competition between different phases in a strongly correlated nano-structure with an 
edge. Making use of the self-consistent Green's function and density matrix renormalization group 
methods, we study a system described by the t-J^ and t-J models on a strip of a square lattice with 
a linear hole density n||. At intermediate interaction strength J/t we find edge stripe-like states, 
reminiscent of the bulk stripes that occur at smaller J/t. We find that stripes attach to edges more 
readily than hole pairs, and that the edge stripes can exhibit a peculiar phase separation. 
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It is well established by now that in the ground state 
of the two-dimensional (2D) strongly correlated cuprates, 
much studied in the context of the high-Tc problemji*^ 
charge carriers tend to be distributed inhomogeneously. 
Such inhomogeneities occur as an attempt to reduce the 
frustration between the exchange and kinetic energies 
arising from the antiferromagnetic (AF) background. ^ It 
is conceivable that such frustration can be further re- 
duced at a boundary. This would induce an effective 
attraction of the bulk inhomogeneity to the edge, thus 
leading to a novel, many-body type edge state. A bulk 
doped Mott insulator might have various types of striped, 
paired, or unpaired ground states. Ordinarily one is 
not concerned with surface states. However, for cuprate 
nanoscale structures, which are of interest due to recent 
progress in nano-synthesisji one may potentially have a 
system with a finite doping in which the holes are bound 
to the edges. In this case, understanding the possible 
edge states is crucial. 

Here, we consider the simplified case of a half-plane, or 
similarly, one end of a long open cylinder, with a vanish- 
ingly small doping sufficient to form a finite linear con- 
centration n|| = Nfi/L of holes near the edge. Nh is the 
number of holes and L is the linear size of the system. 
Given this configuration, several questions arise. For ex- 
ample, if the holes form a single stripe in the bulk, does 
this stripe attach itself to the edge? Do single pairs at- 
tach to the edge? Can one have a one-dimensional (ID) 
"edge phase separation" with the holes clustered together 
in edge droplets? 

In this paper we address some of these questions in 
the context of the t-Jz and t-J models, using the self- 
consistent Green's function and the density matrix renor- 
malization group (DMRG) methods. We obtain an ap- 
proximate edge phase diagram in which the linear hole 
density ny is a control parameter. We find analytically, 
for the t-Jz model, that the stripe edge states do exist in 
a substantial region of the phase diagram. This is verified 
for both the t-Jz and t-J models using DMRG for sys- 
tems up to 11 X 8 sites with periodic/open and 20 x 5 with 
open boundary conditions. We also find that the Ising 
anisotropy enhances tendency to the edge states forma- 
tion. A remarkable feature of these states is that they are 
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FIG. 1: Phase diagram of a system with an edge, in the linear 
hole concentration (ny) v.s. attraction to the edge (J) plane, 
(a) band-like case, (b) case of complex bulk/edge states. The 
labels of the phases indicate where the holes reside. 



reminiscent of and, as we show in some detail, are directly 
related to the bulk stripe states observed in the studies 
of strongly correlated models as well as in real high-Tc 
materialsjiiSii We argue that such edge states can be a 
common feature of strongly correlated nano-systems. 

With our Fig. ^we outline two generic possibilities for 
the formation of the edge state as a function of attraction 
to the boundary (controlled by J) and linear hole con- 
centration Fig. n^a) corresponds to the case when 
both the 2D bulk and the ID edge state are band-like. 
At small J the kinetic energy in the bulk is lower and 
the holes occupy the bulk states only (bulk). At J = Jc 
the bottom of the edge-band becomes degenerate with 
the bulk band. As the attraction to the edge is increased 
further, a fraction of holes from the bulk populates the 
edge band (mixed) . Upon the further increase of J all the 
holes are in the edge band (edge). Note that in the bulk 
and mixed states the chemical potential is pinned at the 
bottom of the bulk band as the finite linear concentra- 
tion n|| = Nfi/L corresponds to zero bulk concentration 
X = Nh/L"^, as L — > cxo. Fig. ^h) outlines a different 
generic scenario in which the edge state formed at J = Jc 
accumulates a finite linear concentration of holes. That 
is, at Jc a finite amount of holes njj becomes attached to 
the boundary at once. This could occur, for example, if 
the lowest energy state is a stripe attached to the edge. 

In this work we focus on a particular case of the 
strongly correlated models, the t-Jz and the t-J models, 
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and, therefore, on a particular case of such a frustration 
characteristic to the doped antiferromagnets. We find, by 
means of analytical and numerical approaches, that the 
edge states do form in these systems and that their phase 
diagram belongs to the class described in Fig. H^b). 

First, we provide a general reasoning for the edge 
states to occur in the i-J-like models. Excitations in 
such models doped with few holes are described as spin- 
polarons, which either form moderately bound pairs or 
move independently.- It is clear that at some J/t such 
excitations would reside close to an edge where they frus- 
trate the AF background less. Thus, edge holes or edge 
pairs can be formed. At larger densities the holes are 
known to be in the "bulk" quasi-lD stripes with internal 
hole concentration close to n\\ — 1/2* Although the ki- 
netic energy advantage for holes in the stripe is substan- 
tial, the magnetic energy cost of the antiphase-domain 
wall (ADW) associated with the stripe is considerable 
too, Ejnag J per holc. However, this magnetic energy 
of the domain wall can be significantly reduced at the 
edge. Therefore, an edge stripe can be formed. One may 
also expect that the low-density stripes will be attracted 
to the edge more readily because their magnetic energy is 
higher {Emag ~ •^/'^||) and its reduction should be more 
effective for them. Thus, one expects that in some range 
of n|| and J an edge stripe can be formed. 

Note that within this scenario edge stripes constitute 
many-body states formed due to a non-trivial collective 
effect. That is, the energy is reduced due to attraction 
of the stripe as an entity. This means that if the ground 
state of the t-J model would be given by a dilute gas 
of spin polarons then the physics outlined in Fig. Qla) 
should be applicable to the bulk and edge bands of these 
excitations. However, if the ground state is a stripe, then 
we expect the edge phase diagram of Fig. Qlb) . 

Another important question here is: what is the 
ground state of the t-J model in the limit of small 
n|| ^ 1/2? The most likely answer is that the ground 
state is given by a closed "stripe loop" with an AF an- 
tiphase shift across the loop and internal hole concentra- 
tion close to the "optimal" n\\ = 1/2, as supported by 
the observation in our previous work, Ref. 0. Such loops 
can be seen as nucleations of the "straight" stripes. In 
this regime we expect the edge states to consist of short 
loop-like segments. 

To demonstrate the described tendencies explicitly and 
quantitatively we consider the t-J model: 

on a strip of a square lattice and in the standard nota- 
tions of constrained fermion and spin operators, J is the 
superexchange and t is the kinetic energy, and summa- 
tion is over nearest-neighbor sites. 

First, we consider a simpler problem where we impose 
a local constraint on the linear hole density. Namely, 
a given n\\ will be assumed to be distributed homoge- 
neously along the boundary. This consideration will help 



to clarify the details of the competition between different 
bulk and edge states such as a dilute gas of holes/pairs 
and "straight" stripes, but will exclude the ID phase sep- 
arated states such as loops. In reality, such a constraint 
may be imposed by the confinement effect of finite sizes 
of a structure or by extra interactions beyond the con- 
sidered models. In a nano-size system single-hole and 
single-pair states can be of interest themselves. 

We will be mostly concerned with the study of the 
anisotropic version of the model, the t-Jz model, for 
which the bulk stripe, pairing, and Nagaoka states, very 
similar to those of the isotropic t-J model, have been 
found in the pastil One can argue that the t-Jz model 
should reproduce qualitatively the holes' ground state 
of the isotropic t-J model, although the details of such 
states might differ. Moreover, for the bulk stripe phase in 
this model a very close quantitative agreement between 
the analytical and numerical methods have been demon- 
strated in our previous work.^ Here we provide an exten- 
sive analysis of the stripes and other excitations of the 
model (QJ in the presence of the boundaries. 

The stripe occurs due to the lowering of the kinetic 
energy of holes at the ID ADW. Although the stripe is 
always beneficial from the point of view of kinetic energy, 
there is a magnetic energy of the ADW, ~ LJz, which 
one needs to compensate. Ignoring hole-hole interaction 
the energy of the bulk stripe is given by: 

^ = :| (n^i - l) + (-11)-^ E, dk , (2) 

where Etotai is relative to Eq of the hole-free Ising AF 
background, fc is a ID wave-vector along the stripe. The 
first term in ^ is the magnetic energy and the second 
term is the energy of the effective ID band of holes filled 
up to the Fermi momentum kp — 7rn||. The quasiparti- 
cle energy Ek is defined through the Dyson's equation: 
Ek -\- 2t cos(fc) — Yi{Ek) — where S(e) is the self-energy 
due to the string- like hole motion away from stripe. 
can be written in a compact form and Dyson equation 
can be solved numerically (see Ref. 0). The bulk stripe 
energy is shown in Fig. Ufa) as a function of n\\ for two 
representative Jz/t values (solid lines). There is an ex- 
cellent agreement of the theory with the DMRG data for 
Jz/t = 0.35 in an 11 X 8 system (circles). 

This bulk stripe energy should be compared to the one 
of bulk holes/pairs. In a homogeneous AF holes attract 
each other in a bound state with the binding energy Ej, 

(~ -Jz/2 for Jz/t > 1 and J^/t"^ for Jz/t < 1).^ 

For the dilute system we assume that the hole pairs do 
not overlap and -Efo/Iz = Nh{Epoiaron + Eb/2), shown by 
the straight lines in Fig. [2Ia). 

In the presence of the boundary two types of new states 
are possible, the edge pair and the edge stripe. Our calcu- 
lations within the retraceable path approximation show 
that the transition from the bulk pair to the edge pair 
occurs at Jz/t w 0.4. For the edge stripe the longitudi- 
nal kinetic energy due to motion along the edge should 
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FIG. 2: (a) Energy of the system per hole as a function of hole 
density n| | for the bulk and edge stripe, and the bulk pair state 
in the t-Jz model (solid, dashed, and dotted, respectively). 
Black circles are the DMRG data for the bulk stripe, Ref. 3. 
Jz/t = 0.35, and Jz/t = 0.1 (inset). Open symbols denote 
transitions: I - pair to edge stripe, II - edge to bulk stripe, 
III - pair to bulk stripe, (b) Constrained (see text) phase 
diagram of the t-Jz system with an edge. A - bulk pairs, B 
- bulk stripes, C - edge stripes, D - edge pairs. 



remain the same as for the bulk stripe while the mag- 
netic and transverse kinetic energy change. At the stripe 
the hole excitations correspond to geometric kinks in the 
magnetic domain wall4i Therefore, placing the stripe at 
the edge eliminates one-half of the broken bonds in the 
ADW and the magnetic energy of the stripe is reduced to 
Jz/2 (l/2n|| — l) per hole. The transverse kinetic energy 
is modified too, which leads to a different form of the 
self-energy T,{e) in the Dyson equation. As a result the 
total energy of the edge stripe is given by the same ex- 
pressions as in Eq. ^ with the reduced magnetic energy 
and different S(£). The resulting Ef^t^i are shown in Fig. 
|2Ia) by the dashed lines. One can clearly see the l/j^n 
behavior in the bulk and edge stripe energies at low n|| 
where the magnetic energy dominates. For Jz/t = 0.35 
the first transition, which occurs at a lower ny « 0.22, is 
from the bulk pairs to the edge stripe state. At higher 
density n|| « 0.58 a transition to the bulk stripe occurs. 
For the smaller value of Jz/t = 0.1 the edge stripe never 
materializes in the ground state and the transition is only 
between the bulk stripe and bulk pairs state (Fig. |2Ia), 
inset). Varying the values of Jz/t in such an analysis wc 
obtain the "constrained" phase diagram in Fig. |2Ib). 

At small Jz/t the boundary is irrelevant and the com- 
petition is between the bulk striped and paired states 
only. At low n|| stripes cannot be formed as Emag oc 
Jz/n\\ is too high. At intermediate Jz/t the edge stripes 
form readily, while it is only at larger Jz/t when the edge 
pairs are formed. Note that the pairs-to-stripes bound- 
ary in Fig. |2Ib) is only weakly J^/t-dependent. At larger 
Jz/t the boundary is between the edge stripes and edge 
pairs. In principle, there can be more subtle varieties of 
the quasi-lD states at the edge, including density-wave 
and others, which might be classified according to the 



Luttinger liquid g-ology:^ Complicating such a classifica- 
tion is the bulk AF background breaking spin symmetry. 

As Jz/t increases the edge state gets strongly attached 
to the boundary and the energy difference between the 
competing phases in this region becomes rather small. In 
that case, the qualitative distinctions between the edge 
pairs and edge stripes, such as the anti-phase shift of the 
AF order parameter across the stripe, diminishes and one 
is not be able to tell what bulk state a given edge state is 
associated with. Thus, one can expect a crossover from 
the edge stripes to the edge pairs at larger Jz/t. 

We have used DMRG to verify numerically the exis- 
tence of these statesi^ We have studied systems Lx x Ly 
up to 11 X 8 sites, with periodic boundary conditions 
(BCs) in the Ly direction and open BCs in direction 
(representing the edges). No external field was applied 
to the open boundariesi^ii Altogether, we have a close 
agreement between our analytical and numerical results. 
All the principal phases, bulk and edge stripe, and bulk 
and edge pair states are identified. Because of the finite 
sizes of the modeled systems only a discrete set of hole 
densities n|| is available, while Jz/t can be varied con- 
tinuously. We have studied the clusters 11 x 4 with one 
hole, 11 X 6 with two holes, and 11 x 8 with two and 
four holes. The number of states kept ranged from 1000 
to 4000, with as many as 20 sweeps, allowing each run 
to converge nicely. An evolution from the bulk to the 
edge stripe was observed in our 11 x 4 (n|| = 0.25) and 
11x6 (n|| = 0.33) clusters. There we start from the bulk 
stripe at small Jz/t in the center of the cluster(i which 
broadens and becomes degenerate with the edge stripe 
upon increase of Jz/t, and then forms the edge state, see 
Fig. Ola). One can find remnants of the antiphase 7r-shift 
due to the ADW for the edge stripe in Fig. EI a), a clear 
signature that this state remains quantum-mcchanically 
connected to the bulk stripe. If, on the other hand, we 
start from the bulk paired state, (11 x 8, two holes) we 
have a transition to an edge state at around Jz/t = 0.4 
(in agreement with the analytical calculation) , which has 
no stripe-like magnetic correlation and we classify it as 
the edge pair. Another distinction is the energy. One can 
compare the energies per hole of the systems of different 
size with the same number of holes in which the edge 
pair and the edge stripe are formed (11 x 6 and 11 x 8 
with two holes). Generally, such per hole energies are 
lower in the edge stripe state, which indicates that: (i) 
the edge stripe is distinct from the edge pair, and that 
(ii) the edge state will be favored in the thermodynamic 
limit. Upon further increasing of Jz/t the energies of the 
edge states become degenerate as shown by such a com- 
parison, which supports our suggestion of the crossover 
from the edge stripe to the edge pair state at larger Jz/t. 

Also, for the concentration ny = 0.25 we observe no 
transition from paired states to stripes as a function of 
Jz/t, in a close agreement with a flat boundary between 
pairs and stripes in our analytical results in Fig. |2Jb). 

Note here that we also found similar states in the case 
of the fully isotropic t-J model, with the phase bound- 
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(a) IJi = 0.45, 2 holes (b) IJt = 0.35, 4 holes (c) IJt = 0.45, 6 holes 



FIG. 3: DMRG results for the edge stripes in the t-J^ model in (a) 11 x 6, (b) 11 x 8, and (c) 20 x 5 clusters. J^/t and the 
number of holes is as indicated. In (a) and (b) only the right half of the cluster containing all the holes is shown with the 
rightmost column being the open boundary. In (c) the bottom row is the "free" open boundary (see text). The arrows represent 
the on-site value of (S?) and the circles are for the hole density (rii). Scales for and (rii) in (b) and (c) are the same. 



aries shifted to somewhat higher J/t and In partic- 
ular, 11 X 8 system with four holes (n|| — 1/2) definitely 
forms an edge stripe at J/t — 0.6. The same system 
exhibits a bulk stripe at J/t = 0.4. For J/t — 0.6, on 
an 11 x 8 cluster, we find the following energies per hole 
for various dopings, in units of t: one hole in the bulk - 
0.90(1), at the edge -0.94(1); one pair, bulk -1.03(1), edge 
-1.03(1); nil = 1/2 stripe, bulk -1.08(1), edge - 1.08(1). 
Despite the degeneracies for the pair and stripe in the 
bulk versus the edge, the edge pair and stripe appear 
quite stable in the DMRG run, and strongly bound to 
the edge. J/t = 0.6 appears to be quite close to the 
transition between bulk and edge stripes and pairs. 

Thus far our DMRG results for the low hole numbers 
have agreed with the "constrained" analytical results in 
which only homogeneous hole states were taken into con- 
sideration. While this analysis can be directly relevant to 
certain nano-scale structures doped with few holes, the 
question remains which phases will survive in the ther- 
modynamic limit. Our energy-per-hole analysis shows 
that energy is always lower in the stripe phase (bulk or 
edge). This implies an instability towards a phase sepa- 
ration of a peculiar kind. That is, the system will tend 
to segregate into a hole empty region and a region where 
a stripe-like phase will be formed. Gas of pairs (edge or 
bulk) will not be selected for the ground state. In fact, 
we observe an indication of such a phase separation al- 
ready in our t-Jz 11x8 cluster with four holes (ri|| = 0.5). 
We start at a lower J^/t with the "regular" bulk stripe 
which evolves into an excellent edge stripe at larger Jz/t^ 
similar to one in Fig. Eta). However, on its way to the 
edge it forms a loop-like stripe object, with clear ADW, 
attached to the edge with its ends, see Fig. ^h). 



We augmented our study of higher hole concentrations 
with the 20 X 5 cluster with open BCs for all the edges 
and a strong repulsive potential for holes (in the form of 
a chemical potential fi = +t) along the left, upper, and 
right edges to create a preferential boundary and look 
for potential edge states. Hole numbers up to six were 
studied. They organized themselves into an edge arc-like 
stripe, see Fig. OIc). Similar patterns have been pro- 
duced by four holes as well. Such states show clear sig- 
natures of the ADW and are very stable. An "internal" 
hole concentration per arc unit length can be estimated 
and is close to n|| = 1/2. Thus, one can argue that the 
ground state of the semi-infinite t — J system is likely to 
be composed of either "straight" or "loop" stripes (de- 
pending on a given n||), which attach themselves to the 
boundary at some J/t. Note, that nowhere in our study a 
mixed bulk-edge state was detected and, as we discussed, 
the t-Jz or t-J system with edges should belong to the 
class described in our Fig. ^h). 

In summary, we have studied the edge phase dia- 
gram of a strongly-correlated nano-structure and found 
that it consists of the bulk pairs, bulk stripes, and edge 
stripes. Using DMRG and analytical methods we have 
shown that these edge stripe states are directly related 
to the bulk stripe states and have their origin in the non- 
trivial interplay between magnetic and kinetic energies 
in strongly-correlated systems. 
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